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RESUMEN 

El programa GID es un pre y post procesador que facilita la entrada y salida de 
resultados, en forma grafica, de estructuras analizadas con elementos finitos. En este artfculo 
se presenta el uso de este programa como post procesador aplicado al analisis sfsmico de la 
Presa la Esperanza que se halla ubicada en la Provincia de Manabf. 

Se indica como se generan los archivos Flavia.msh y Flavia.res. En el primer archivo 
se genera la malla de elementos finitos y en el segundo todos los resultados que se desean 
graficar. Este artfculo es un complemento al uso del programa GID que fue publicado por los 
autores en el segundo volumen de la Revista Cienciade 2010. 

ABSTRACT 

The program denominated GID is a pre and post processor that makes the input and 
output of fata easier, of structures analyzed graphically with finite elements. This article presents 
the use of this program as a post processor applying the seismic analysis of La Esperanza Dam, 
which is located in Manabf Province. 

It is indicated how to generated the Flavia.msh and Flavia.res files. The finite element 
mesh is generated in the first file and all the results are generated in the second file, wich are 
wanted to be plotted. This article is a complement for the use of the program GID that was 
published by its authors in the second edition of Science Magazine in 2010. 
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1. INTRODUCCION 

En Aguiar et al. (2010) se indico la forma de usar el Programa GID como pre y post 
procesador pero se dio mas entasis a la primera opcion. Ahora se va a explicar con mas detalle 
la segunda opcion considerando que el usuario del GID va a resolver analizar cualquier 
problema con programas de CEINCI-LAB de elementos finitos y quiere presentar los 
resultados en forma elegante, agradable y comprensiva, que se logra con el Programa GID. 

Para cumplir con este proposito se deben generar los archivos para el Post-proceso, 
archivo flavia.msh y flavia.res a traves del MATLAB. 

Para un mejor entendimiento se realiza un ejemplo de una presa de hormigon a 
gravedad, la misma que se muestra en la figura 1 la parte de la izquierda, con color azul 
corresponde al agua y la parte de la derecha a la Prera cuyas dimensiones se indican en metros 
en la figura 2. 



Figura 1 Presa de hormigon a gravedad. 


En primer lugar se realiza el analisis sfsmico de la Presa empleando el Metodo de 
Superposition Modal utilizando un espectro de respuesta inelastico. Aguiar (2008). Luego se 
encontrara la Presion Hidrodinamica para el coeficiente sfsmico que se obtiene del analisis 
sfsmico. 

Para el analisis sfsmico de la Presa se trabaja con el elemento finito Q4, que considera 
dos grados de libertad por nodo (desplazamiento horizontal y vertical) pero este elemento finito 
es suavizado por el efecto de flexion. Onate (1992). Se consideran dos divisiones en sentido X y 
dos divisiones en sentido Y, de tal manera que se trabaja con una malla de 4 elementos finitos. 
Se consideran tan pocos elementos finitos por facilidad de la explication. 

La base de la Presa se encuentra restringida a los grados de libertad (hipotesis de 
calculo) por lo que el sistema posee 6 grados de libertad en la direccion X y 6 en la direccion Y, 
los mismos que se presentan a la derecha de la figura 3. Notese que primero se han numerado 
de corrido todos los grados de libertad en sentido horizontal y posteriormente todos los grados 
de libertad en sentido vertical esto es con el proposito de poder trabajar con la matriz de rigidez 
condensada en el analisis sfsmico. Aguiar (2004, 2007). A la izquierda de la figura 3 se presenta 
la numeration de los nodos y de los elementos. 



Presentation grafica en el GID a partir de resultados obtenidos en MATLAB 


155 



Figura 2 Dimensiones de la Presa en metros. 



Figura 3 a) Elementos de la malla y nodos. b) Grados de libertad. 


2. ARCHIVO FLAVIA.MSH 

Contiene los datos para generar la geometrfa de la malla, para lo cual requiere que se 
indiquen las coordenadas de cada nodo ( coordenadas ) y la conectividad de cada elemento 
finito (conectividad). En parentesis se ha indicado el nombre de los programas de la librerfa 
CEINCI-LAB que se explican a continuation. 

> Uso de Programa coordenadas 

Este programa encuentra las coordenadas de los nodos de la malla de elementos finitos. 
Para la estructura de la figura 2, con la numeracion de los nodos de la figura 3a. Estas 
coordenadas son las indicadas en la tabla 1 . La forma de uso del programa es: 

[ coord]=coordenadas(altura, divx, divy,alfa, beta, base) 

Donde altura es la altura de la Presa; divx es el numero de divisiones en sentido X; 
divy es el numero de divisiones en sentido Y que se consideran para generar la malla; alfa es 
la relacion de la longitud aguas arriba de la Presa con la altura; beta es la relacion de la longitud 
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aguas abajo de la Presa, con la altura; base es la dimension inferior de la Presa. El programa 
en la variable coord reporta las coordenadas en sentido X, Y de los nodos de la Presa. Para la 
Presa analizada se tiene: [coord]=coordenadas(11,2,2,3/11,7/11,13) 


Tabla 1 Coordenadas de los nodos. 


Nodo 

X 

Y 

1 

0 

0 

2 

6.5 

0 

3 

13 

0 

4 

1.5 

5.5 

5 

5.5 

5.5 

6 

9.5 

5.5 

7 

3 

11 

8 

4.5 

11 

9 

6 

11 


function [coord] ^coordenadas (altura, divx, divy, alf a, beta, base) 

longy=altura/divy; 

for i=0 : divy 

Bases (i+1) = (base- (altura/divy) *i*alfa- (altura/divy) *i*beta) ; 
xvec (i+l)=(altura/ divy) *i*alf a; 

end 

for i=l : divy+1 

xvecm ( i ) =xvec ( i ) +Bases ( i ) / 2 ; 
xvecd ( i ) =xvecm ( i ) +Bases ( i ) / 2 ; 

end 

for i=l : divy+1 

for j = 1 : divx+1 
aux=l ; 

cxl (i, j ) =xvec (aux*i) + (Bases (aux*i) / divx) * ( j -1) ; 
aux=aux+l ; 
end 

end 

for i=l : divy+1 

for j = 1 : divx+1 

cyl (i, j ) =longy* (i-1) ; 

end 

end 
aux=l ; 

for i=l : divy+1 

for j =1 : divx+1 
cx (aux) =cxl (i, j ) ; 
cy (aux) =cyl (i, j ) ; 
aux=aux+l ; 
end 

end 

coord= [ cx (:,:); cy (:,:)]' ; 

% fin 
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> Uso de Programa conectividad 

Este programa muestra la conectividad de cada elemento finito en forma antihoraria, 
empezando en el nodo inferior izquierdo. Se usd la programacion de la subrutina NUDE de 
CEINCI-LAB, que para este estudio se la denomino conectividad. El programa genera la 
matriz denominada CONEC que contiene la conectividad de los nodos. En la tabla 2 se indica 
esta matriz y la forma de uso es la siguiente: [conec]=conectividad(2,2) 


Tabla 2 Conectividad de los nodos. 


Elemento finito 

Nodo i 

Nodo ii 

Nodo iii 

Nodo iv 

1 

1 

2 

5 

4 

2 

2 

3 

6 

5 

3 

4 

5 

8 

7 

4 

5 

6 

9 

8 


function [conec] ^conectividad (divx, divy) 
nef i=divx*divy ; 

% NUDOS DE CADA ELEMENTO FINITO 
aux2=l ; 
aux3=l ; 
for i=l : divx 

NUDE ( i , : ) = [ aux2 , aux2+l, aux2+l+divx+l , aux2+l+divx] ; 
aux2=aux2+l ; 
end 

for j = 1 : divy-1 

for i= (divx+1 ) * j : nef i 

NUDE (i, : ) =NUDE (i-divx, : ) + (divx+1) ; 

end 

end 

conec=NUDE ; 

% fin 


> Uso de Programa GID_mesh 


Con las matrices que contiene a las: coordenadas (coord.) y la conectividad (conec) se 
usa el programa GID_MESH (Mrogisnki, 2010), el mismo que genera de manera automatica el 
archivo flavia.msh. A continuacion se lo presenta con unas modificaciones realizadas para el 
presente estudio: 


function GID_mesh (problem, coord, conec, divx, divy) 

whole_problem = [ problem '.flavia.msh']; 
fid = f open (whole_problem, ' w ' ) ; 


npoin= (divy+1 ) * (divx+1 ) ; % Numero de puntos o nodos de la malla 
nelem=divy*divx; % Numero de elementos finitos 


linea ( 1 , : ) 
linea (2 , : ) 
linea (3, : ) 
linea ( 4 , : ) 
linea (5, : ) 


MESH dimension 2 ElemType Quadrilateral Nnode 4 ' ; 
Coordinates ' ; 
end coordinates ' ; 
Elements ' ; 
end elements ' ; 


%ESCRIBE ENCABEZADO: PRIM. LINEA 
fprintf (fid, '%s\n',linea(l, :) ) ; 
%COORDENADAS 
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fprintf (fid, ' %s\n ' , linea (2 , : ) ) ; 
coordenadas= [ (1 : npoin) ' coord] ; 

fprintf (fid, ' %5 . Of %19.7f %1 9 . 7f \n ' , coordenadas ' ) ; 
fprintf (fid, , %s\n , ,linea(3, :) ) ; 

%CONECTIVIDADES 

fprintf (fid, , %s\n , ,linea(4, :) ) ; 
conectiv= [ (linelem) ' conec] ; 

fprintf (fid, 1 %5 . Of %5.0f %5.0f %5.0f %5 . Of \n ' , conectiv ' ) ; 
fprintf (fid, , %s\n , ,linea(5, :) ) ; 

f close (fid) ; 

return 


La forma de uso del programa GID_mesh es la siguiente: 

GID_mesh(problem, coord, conec, di vx, divy) 

Donde problem es el nombre del programa; coord es el nombre de la matriz que 
contiene a las coordenadas en sentido X, Y; conec es el nombre de la matriz con la 
conectividad de los elementos finitos; divx es el numero de divisiones en sentido X; divye s el 
numero de divisiones en sentido Y. Para este ejemplo a la Presa en estudio se la ha 
denominado ‘presa_hormigon’. Por lo tanto, presa_hormigon.flavia.msh es el archivo que se 
generara el mismo que contiene la siguiente information: 

MESH dimension 2 ElemType Quadrilateral Nnode 4 


Coordinates 

1 

0.0000000 


0.0000000 

2 

6.5000000 


0.0000000 

3 

13.0000000 


0.0000000 

4 

1.5000000 


5.5000000 

5 

5.5000000 


5.5000000 

6 

9.5000000 


5.5000000 

7 

3.0000000 


11.0000000 

8 

4.5000000 


11.0000000 

9 

6.0000000 


11.0000000 

end coordinates 



Elements 
1 1 

2 5 

4 


2 2 

3 6 

5 


3 4 

5 8 

7 


4 5 

6 9 

8 


end elements 




3. ARCHIVO FLAVIA.RES 

Este archivo contiene los resultados obtenidos con CEINCI-LAB. Pero para ser 
visualizados de la forma deseada se deben tomar en cuenta los siguientes aspectos. 

• El GID muestra los resultados por nodo. Por lo tanto se debe ubicar los grados de 
libertad en el nodo respectivo. (Ver tabla 4). 

• Los nodos restringidos (1 , 2 y 3) se deben rellenar con cero como se hizo en la tabla 4. 

Debido a estas consideraciones de deben modificar, muy brevemente los resultados que 
se obtienen con los programas de CEINCI-LAB para adaptarlos al programa GID_result para 
que este acorde con los valores de la tabla 4. 
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Para el ejemplo, las propiedades del hormigon de la estructura analizada son las 
siguientes: Modulo de Elasticidad 2100000 T/m 2 ; coeficiente de Poisson v -0.2. El sismo de 
analisis es el espectro de diseno inelastico hallado para la Presa: Minas-La Union- Pitahuina, 
para el sismo MCE (Maximo Sismo Crefble, por sus siglas en ingles) que tiene un perfodo de 
retorno de 970 anos. Aguiar (2010). 

Como se indico anteriormente, se aplico el Metodo de superposition modal aplicando el 
criterio de combination modal de la Normativa Sfsmica de Peru 2003 y se obtuvieron los 
resultados indicados en la tabla 3. 

Tabla 3 Desplazamientos inelasticos y Fuerzas sismicas 


GDL 

qine [m] 

Ft [T] 

1 

1.15E-04 

17.088 

2 

1.09E-04 

33.869 

3 

1.08E-04 

14.666 

4 

3.43E-04 

19.074 

5 

3.42E-04 

37.874 

6 

3.22E-04 

17.403 


En la tabla 3; qine es el vector de desplazamientos inelasticos hallados a partir de los 
desplazamientos elasticos multiplicados por el factor de reduccion de respuesta de las fuerzas 
sismicas R . Para este ejemplo se trabaja con R - 1 pero para el siguiente ejemplo que 
corresponde a la Presa La Esperanza se trabaja con R = 2 . La otra variable de la tabla 3 es Ft 
que contiene a las fuerzas maxi mas probables debido a sismo. 

Como se observa, CEINCI-LAB da los resultados en cada grado de libertad. Para 
poder graficar en el GID estos valores es necesario reordenarlos tomando como referencia el 
nodo, como se ilustra en la tabla 4. 


Tabla 4 Desplazamientos y Fuerzas para e 


Nodo 

qine gid 
[m] 

Ft_gid [T] 

1 

0 

0 

2 

0 

0 

3 

0 

0 

4 

1.15E-04 

17.088 

5 

1.09E-04 

33.869 

6 

1.08E-04 

14.666 

7 

3.43E-04 

19.074 

8 

3.42E-04 

37.874 

9 

3.22E-04 

17.403 


GID 


Es importante destacar que los desplazamientos en los nodos 1 , 2 y 3; efectivamente 
son cero debido al modelo indicado en las figuras 2 y 3 (Apoyo Fijo). Pero las fuerzas en los 
nodos 1, 2 y 3, no son cero ya que son las reacciones horizontales que se generan en el 
apoyo. Pero el objetivo del artfculo esta orientado al uso grafico del GID por este motivo se 
pasa por alto esta consideration. 

Los programas de CEINCI-LAB generan los archivos qine y Ft. Pero estos archivos no 
tienen el valor en todos los nodos (falta de los nodos restringidos). Por este motivo se debe 
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crear unos nuevos archivos con la terminacion “_gid”. Se puede colocar cualquier otra 
terminacion en el nombre. Con esta acotacion y para estar en concordancia con lo indicado en 
la tabla 4 se debe proceder de la siguiente manera: 

% Desplazamientos modales inelasticos 
qine_gid= [ zeros (divx+1 , 1 ) ; qine] ; 


% Fuerzas Sismicas 
Ft_gid= [ zeros (divx+1 , 1) ; Ft] ; 


Siguiendo con el analisis, se hallaron las fuerzas hidrostaticas e hidrodinamicas, a las 
que se anaden las fuerzas sismicas y el peso propio, encontrando un vector de cargas totales. 
Con este vector de cargas se calcula los desplazamientos totales que se muestran en la tabla 5. 
Las fuerzas hidrodinamicas se encontro utilizando el modelo de Housner (1978) descrito en 
Chwang y Housner (1978). 


T abla 5 Desplazamientos totales 


GDL 

q [m] 

1 

1.67E-04 

2 

1.59E-04 

3 

1.52E-04 

4 

3.97E-04 

5 

3.99E-04 

6 

3.85E-04 

7 

4.09E-05 

8 

-6.74E-05 

9 

-1.08E-04 

10 

1.13E-05 

11 

-6.70E-05 

12 

-1.32E-04 


El vector q contiene los desplazamientos horizontales y verticales en los nodos de la 
Presa, pero para poder graficar hay que separar cuales son los desplazamientos en sentido X y 
cuales son los desplazamientos en sentido Y. Pero a mas de ello se debe colocar ceros en los 
nodos restringidos esto se lo hace de la siguiente manera. 

% Desplazamientos totales 
q_gid= [q ( 1 : na, 1 ) , q (na+1 : ngl , 1 ) ] ; 
q_gid= [ zeros (divx+1 , 2 ) ; q_gid] ; 

% na=grados de libertad principales, en este caso son 6 

Con ello se obtuvo los resultados de la tabla 6. Puede verse la figura 3 para apreciar los 
nodos y los grados de libertad del sistema. 

Notese que los desplazamientos son colocados en una misma variable, en la columna 
1 los desplazamientos en X y en la columna 2 los desplazamientos en Y. 

Este mismo proceso se realiza para las fuerzas hidrostaticas, hidrodinamicas y 
cualquier resultado que desee graficarse. 
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Tabla 6 Desplazamientos totales para el GID. 


Nodo 

q_g 

d[m] 

X 

Y 

1 

0 

0 

2 

0 

0 

3 

0 

0 

4 

1.67E-04 

4.09E-05 

5 

1.59E-04 

-6.74E-05 

6 

1.52E-04 

-1.08E-04 

7 

3.97E-04 

1.13E-05 

8 

3.99E-04 

-6.70E-05 

9 

3.85E-04 

-1.32E-04 


> Uso de Programa GID_result 

El programa GID_result (Mroginski, 2010), es el encargado de generar 
automaticamente el archivo flavia.res el mismo que contiene todos los resultados que se 
desean graficar. Para lo cual llama a todas las nuevas variables (qine_gid, Ft_gid, q_gid,...)- A 
continuation se presenta dicho programa el cual ha sido adecuado para el presente estudio: 

function 


GID result (problem, q gid,Q gid, FH gid, FHI gid. Ft gid, qine gid, 

y) 

whole problem = [ problem ’.flavia.res']; 
fid = f open (whole problem, ' w ') ; 

Peso 

_gid. 

divx, div 

o 

npoin= (divy+1) * (divx+1) ; 





o 

linea ( 1 , 

) = [’DESPLAZAMIENTOS 1 

1 

2 

1 

O' ] ; 

linea (2 , 

) = ['VECTOR CARGAS 1 

1 

2 

1 

0 ' ] ; 

linea (3, 

) = [ 'HIDRODINAMICAS 1 

1 

2 

1 

0']; 

linea ( 4 , 

) = [ 'HIDROSTATICAS 1 

1 

2 

1 

0']; 

linea (5, 

) = [ ’ FUERZA_SISMICAS 1 

1 

1 

1 

O' ] ; 

linea ( 6 , 

) = ['DESP INELASTICO 1 

1 

1 

1 

0']; 

linea ( 7 , 

) = [ ’ PESO_PROPIO 1 

1 

1 

1 

O' ] ; 

% La Tercera columna indica el numero de columnas del 

archivo 




% DESPLAZAMIENTOS TOTALES 





fprintf (fid 

, '%s\n’,linea(l, :) ) ; 





desptotal= [ 

( 1 : npoin) ' q gid] ; 





fprintf (fid 

CL 

, '%5.0f %23.15f %23 . 15f \n ' , desptotal ' ) ; 





O 

% VECTOR CARGAS 





fprintf (fid 

, '%s\n’,linea(2, :) ) ; 





fuerzas= [ (l:npoin)' Q gid]; 





fprintf (fid 

, '%5.0f %23.15f %23.15f\n',fuerzas') ; 





"o 

%FUERZAS 

HIDRODINAMICAS 





fprintf (fid 

, '%s\n',linea(3, :) ) ; 





dinamica= [ 

(1: npoin)' FH gid]; 





fprintf (fid 

CL 

, '%5.0f %23.15f %23 . 15f \n ' , dinamica ' ) ; 





O 

%FUERZAS 

HIDROSTATICAS 





fprintf (fid 

, '%s\n',linea(4, :) ) ; 





estatica= [ 

(1: npoin) ' FHI gid] ; 





fprintf (fid 

, '%5.0f %23 . 15f %23.15f\n', estatica ' ) ; 
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%FUERZAS SISMICAS 
fprintf (fid, ' %s\n ' , linea (5, : ) ) ; 
sismica=[ (l:npoin) ' Ft_gid] ; 
fprintf (fid, ’ %5 . Of %23.15f\n f , sismica ’ ) ; 

" 6 

%DESPLAZAMIENTOS MODALES INELASTICOS 
fprintf (fid, , %s\n , ,linea(6, :) ) ; 
despine= [ (l:npoin) ' qine_gid] ; 
fprintf (fid, '%5.0f %23.15f\n', de spine ’ ) ; 

" 6 

%PESO_PROPIO 

fprintf (fid, '%s\n',linea(7, :) ) ; 

peso=[ (l:npoin)' PESO]; 

fprintf (fid, * %5 . Of %23.15f\n' ,peso' ) ; 


f close (fid) ; 
return 


Como puede darse cuenta el programa llama a nuevas variables, que son las siguientes: 

•/ q_gid = Desplazamientos totales. 

Q_gid = Vector de Cargas Totales. 

V FH_gid = Fuerzas Hidrodinamicas. 

S FHI_gid = Fuerzas Hidrostaticas. 

^ Ft_gid = Fuerzas Sfsmicas. 

S qine_gid = Desplazamientos inelasticos. 

S Peso_gid = Peso propio. 


Ahora se explica un aspecto importante del programa GID_result que hay que tomar en cuenta 
en dos partes del mismo: 


• Primera parte. 


linea ( 1 , 

) = ['DESPLAZAMIENTOS 

1 

1 

2 

1 

O' ] ; 

linea (2 , 

) = ['VECTOR CARGAS 

1 

1 

2 

1 

0']; 

linea (3, 

) = ['HIDRODINAMICAS 

1 

1 

2 

1 

0']; 

linea ( 4 , 

) = ['HIDROSTATICAS 

1 

1 

2 

1 

0 ' ] ; 

linea (5, 

) = [ FUERZA_SISMICAS 

1 

1 

1 

1 

O' ] ; 

linea ( 6 , 

) = [ ' DESP_INELASTICO 

1 

1 

1 

1 

0 ' ] ; 

linea ( 7 , 

) = ['PESO PROPIO 

1 

1 

1 

1 

0']; 


Puede observarse que existen 5 numeros a la derecha en cada linea, el tercer termino en 
unos casos vale 1 y en otros vale 2. La razon es porque cuando una variable tiene una 
columna (resultados solo en X) hay que ubicar 1 , y cuando posee dos columnas (resultados 
en X e Y) se coloca 2. 

• Segunda parte. 

% DESPLAZAMIENTOS TOTALES 
fprintf (fid, , %s\n , ,linea(l, :) ) ; 
desptotal= [ (l:npoin) ’ q gid] ; 

fprintf (fid, %23.15f %23 . 15f \n 1 , desptotal 1 ) ; 

o, 

o 

% DESPLAZAMIENTOS MODALES INELASTICOS 
fprintf (fid, , %s\n , ,linea(6, :) ) ; 
despine= [ (linpoin) ’ qine_gid] ; 
fprintf (fid, ' %5 . Of %23 . 15fl\n ', de spine ' ) ; 
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Los digitos resaltados en verde %5.0f son para mostrar los nodos de la malla. No tiene 
decimales. Las cifras resaltadas en amarillo indican el numero de columnas que posee esa 
variable. En el caso de desplazamientos inelasticos, fuerzas sfsmicas y peso propio se coloca 
%23.15f indicando que exista un campo de 23 columnas que contiene 15 decimales (formato) 
de resultados. Pero en el caso de desplazamientos totales, fuerzas hidrodinamicas, 
Hidrostaticas y Cargas totales se coloca %23.15f %23.15f indicando que contiene dos columnas 
de resultados. El archivo que se genera para este ejemplo es presa_hormigon.flavia.res: 


DESPLAZAMIENTOS 

1 0.000000000000000 

2 0.000000000000000 

3 0.000000000000000 

4 0 . 000167370046953 

5 0.000159039820593 

6 0 . 000152194758854 

7 0.000396905995904 

8 0.000398572202643 

9 0.000384654309953 
VECTOR_CARGAS 

1 0.000000000000000 

2 0.000000000000000 

3 0.000000000000000 

4 57 . 041223579405852 

5 53.845670931776596 

6 14 . 666441152647547 

7 19.073700322549364 

8 37 . 873684126802459 

9 17 .403259809418127 


HIDRODINAMICAS 

1 0.000000000000000 

2 0.000000000000000 

3 0.000000000000000 

4 13.369941378536376 

5 6.684970689268188 

6 0.000000000000000 

7 0.000000000000000 

8 0.000000000000000 

9 0.000000000000000 

HIDROSTATICAS 

1 0.000000000000000 

2 0.000000000000000 

3 0.000000000000000 

4 26.583333333333332 

5 13.291666666666666 

6 0.000000000000000 

7 0.000000000000000 

8 0.000000000000000 

9 0.000000000000000 


FUERZA_SISMICAS 

1 0.000000000000000 

2 0.000000000000000 

3 0.000000000000000 

4 17 . 087948867536142 

5 33.869033575841740 

6 14 . 666441152647547 

7 19.073700322549364 

8 37 . 873684126802459 

9 17 .403259809418127 
DESP INELASTICO 


0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000040906785858 
- 0.000067431516141 
- 0.000108408834078 
0 . 000011270915070 
- 0 . 000067021911310 
- 0 . 000131808241413 
112 1 
0.000000000000000 
0.000000000000000 
0.000000000000000 
- 50 . 994521473037608 
- 69.299999999999997 
— 34 . 649999999999999 
- 18.149999999999999 
- 36.299999999999997 
- 18 . 149999999999999 
112 1 

0.000000000000000 
0.000000000000000 
0.000000000000000 
- 5.469521473037609 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
112 1 

0.000000000000000 
0.000000000000000 
0.000000000000000 
- 10.874999999999998 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
0.000000000000000 
1111 


0 


0 


0 


0 


0 


0 
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1 0.000000000000000 

2 0.000000000000000 

3 0.000000000000000 

4 0 . 000114913971891 

5 0.000108500381712 

6 0 . 000107787590213 

7 0 . 000342863423571 

8 0 . 000342251655241 

9 0.000322338015231 


PESO_PROPIO 

1 0.000000000000000 

2 0.000000000000000 

3 0.000000000000000 

4 -34 . 649999999999999 

5 -69.299999999999997 

6 -34 . 649999999999999 

7 -18.149999999999999 

8 -36.299999999999997 

9 -18.149999999999999 


0 


4. GRAFICOS OBTENIDOS CON EL GID 


Una vez que se ejecuta el programa GID, se debe cambiar al Postproceso 
seleccionando el noveno icono que aparece con luz encendida, ver figura 4. 


Archive Ver Geometria Utilidades Dates Malla Calcular Ayuda 


^ LayerO T | ^ ? | 

$ • 

Cambiar entre prey postproceso 

7 



Figura 4 Seleccion del icono para que trabaje como postproceso 


Luego se abre la carpeta (con el segundo icono) en la cual se tienen los programas (se 
busca con la flecha indicadora del cuadro inferior), en este caso se encuentran en la carpeta 
denominada Presa hormigon GID (se selecciona esta carpeta) y se tiene lo mostrado en la 
figura 5. 

Archivo Ver Utilidades Cortar Ver resultados Opciones Ventana Ayuda 









'W A 


& 

/\ 


-V* 

» ^ » 




r 

Leer Postproceso 


Directorio: Presa hormigon GID 

jlI e # 

B presa_hormigon. flavia. msh 


B presa_hormigon. flavia. res 


3 

J 

N ombre del archivo: 

Abrir 

Archivos de tipo: postproceso de GiD (*.res/.msh/.bin) 

T Cancelar 


Figura 5 Seleccion de la carpeta con archivos Flavia. 
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Ahora se selecciona cualquiera de los dos archivos y al seleccionar el icono de ver 
resultados se tiene lo indicado a continuacion y el lector puede empezar a mostrar lo que desea 
ver. En la figura 6 se indica como aparece actualmente la pantalla del GID. 



Figura 6 Despliegue de archivos creados que se desean visualizar. 

Se deja al lector que vea los siguientes resultados, si tiene dificultad puede remitirse al 
artfculo de Aguiar et al. (2010) donde se indica la forma de usar el GID como post proceso. Con 
estas indicaciones en la figura 7 se muestra la numeration de los nodos y elementos de la malla 
de elementos finitos. 



Figura 7 Malla de elementos finitos. 
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En la figura 8 se muestran los desplazamientos laterales ante la accion sfsmica 
considerando R = 1 es decir se tratan de desplazamientos elasticos pero como se puede 
trabajar con un valor de R ^ lse indica que son desplazamientos inelasticos. 


y 



Areas coloreadas de DESP INELASTICO. 


DESP INELAS TICQ 

1 0.00034286 
0.00030477 
0.00026667 
0.00022858 
■ 0.0001 9048 


I 


0.00015238 

0.00011429 

7.6192e-05 

3.8096e-05 

0 


= igura 8 Desplazamientos inelasticos, hallados con R=1 estrictamente son desplazamientos elasticos. 


y 



X-DESPLAZAMIENTC 


Areas coloreadas de DESPLAZAMIENTOS, X- DE SPLAZ AM I ENTOS . 


Figura 9 Desplazamientos totales en X, considerando sismo, presiones y peso propio. 
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En la figura 9 se presentan los desplazamientos horizontales considerando: accion 
sismica, presion hidrostatica, presion hidrodinamica y peso propio. 

En la figura 10 se indica las fuerzas estaticas equivalentes, halladas en cada uno de los 
nodos por efecto de la accion sismica. 



FUERZA SISMICAS 


I 


I 


37.874 

33.665 

29.457 

25.249 

21.041 

16.833 

12.625 

8.4164 

4.2082 

0 


Figura 10 Fuerzas sfsmicas. 


5. APLICACION A LA PRESA LA ESPERANZA 


En Aguiar et al (2010), se realizo el analisis estatico de una columna de suelo confinado 
de La Presa La Esperanza, luego ahf se puede complementar la informacion que se indica en 
este artfculo. En la figura 1 1 se muestra la seccion transversal con sus diferentes materiales. 



Figura 11 Seccion transversal y materiales. 


La Atagufa fue la primera obra que construyeron para desviar las aguas del rio Carrizal, 
fue construida con una arena arcillosa compactada al 100%. Para proteger la Atagufa se coloco 
material enrocado delante de ella, como se aprecia a la izquierda de la figura 1 1 , son bolones 
de rocas con un tamano maximo de 60 cm., Luego se construyo el Nucleo con el mismo 
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material de la Atagufa y finalmente el espaldon que es una grava arenosa. En la figura 2 y en la 
tabla 1 de Aguiar etal. (2010) se presenta con detalle los materiales de la Presa. 

5.1 Modelo Utilizado 

Para realizar el analisis sfsmico se genera una seccion equivalente que represente la 
estructura real de la presa. El objetivo de realizar este proceso es obtener una geometrfa mas 
sencilla para efectuar dicho analisis, ya que la estructura real posee distintos taludes y bermas 
que dificultarfan el calculo. 

En la figuras 12 se presenta la geometrfa real de la Presa y en la figura 13 el modelo 
equivalente utilizado para el analisis sfsmico. 



Figura 12 Seccion real de la Presa La Esperanza. 



Figura 13 Seccion equivalente utilizada en el estudio. 


5.2 Estudio de Suelos 

Se realizaron varios ensayos de suelos, para encontrar las propiedades de los 
materiales que conforman la Presa, uno de ellos fue el que se presenta en la figura 14, en que 
se observa una muestra de arena arcillosa sometida a un ensayo de compresion para 
determinar el modulo de elasticidad del suelo. 

En la parte superior izquierda se muestra al material siendo compactado con el Proctor, 
la muestra se dividio en tres capas y en cada una de ellas se compacto con 25 golpes; la 
humedad natural de la muestra fue de 40%, razon por la cual se dejo secar hasta que tenga una 
humedad del 27% que es la que debe tener el nucleo de la Presa. 

En la parte superior derecha se muestra el momento en que se pesaba la muestra para 
obtener el peso especffico, de paso se observa el material antes de ser ensayado. En la parte 
inferior izquierda se ve la muestra durante el ensayo a compresion; notese que se adapto un 
deformfmetro para encontrar la deformacion vertical y en la parte inferior derecha se aprecia la 
muestra de material ya colapsado, el criterio de colapso se establecio cuando sin un aumento 
de carga la muestra se deforma apreciablemente. 
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Figura 14 Ensayos realizados a una muestra de arena arcillosa. 

En el estudio se vio que el modulo de elasticidad de la arena arcillosa se ajusta a la 
siguiente ecuacion. 


Es=^fs (D 

Donde f s es la resistencia a la compresion de la arena arcillosa en kg/ cm 2 ; E s es 

el modulo de elasticidad del suelo o modulo de Young en kg/ cm 2 . Tambien se obtuvo el peso 
especifico y el coeficiente del Poisson del material que conforma el Nucleo y Espaldon de la 
Presa, estos resultados se indican en la tabla 7. 


Tabla 7 Propiedades de los materiales 


Material 

Peso especifico 
[T/m3] 

Modulo de Young 
[T/m2] 

Coeficiente 
de poisson 

Nucleo 

1.90 

3000 

0.48 

Espaldon 

1.90 

24440 

0.36 


El material de la ataguia no se lo muestra en la tabla anterior porque tiene las mismas 
caracteristicas que el nucleo. 
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5.3 Malla de elementos finitos 

Para obtener la malla de elementos finitos se trabajo con 10 divisiones en X e Y. En la 
figura 15 y 16 se muestran los elementos finitos y nodos de la malla respectivamente. 
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Figura 15 Numeracion de los elementos Finitos. 
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Figura 16 Numeracion de los nodos. 


La estructura en analisis posee 110 grados de libertad en la direccion X y 110 en la 
direccion Y. En la figura 17 se presenta la numeracion de los grados de libertad principales, 
para el analisis sfsmico ante la componente horizontal del suelo. 


100 ^ 101 ^ 10 ^ 10 ^ 10 ^ 10 ^ 106 ^ 107 ^ 108 ^ 109 ^ 110 
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89 90 91 92.93,94,95,! 
78.79 .80 ,81 ,82 ,83,, 84 , 

36 97 98 99 

,85 ,86 ,87 , 88 

67, , 68 ,, 69 ( ,70 ( [ 7l‘ (t 72 , 73 , 

, 74,, 75 ,, 76,, ,77, 

56,, 57, 58, 59 , 60 , 61 , 62 , 

63 , 64 , 65 , 66, 

45,, 46,, 47,, 48 ,, 49 , , 50 , , 51 , 

, 52,, 53,, 54,, 55 f 

34 i. 35j. 36 ,, 37,, 38,, 39,, 40, 

, 41 , , 42 , , 43 , , x 44 ( 

23, , 24 , , A 25 , , A 26 , , 27 ,, 28 , , 29 , 

f 30,. 31 M 32,. * ,33, 

12,, 13,, 14 15 ,, 16 ,, 17 ,, 18_ 

19 ,, 20,, 21,, ,22, 

1 2 ,, 3 ,, 4,, 5,, 6,, 7 

n 8 i. 9 i. * 1 0 , , ,11, 




Figura 17 Grados de libertad principales. 


Se utilizo el espectro inelastico hallado para la Presa Minas-La Union-Pitahuina, 
considerando un valor de R = 2 . 


El objetivo del artfculo es ilustrar la forma como se trabaja el programa GID como 
postprocesador a partir de resultados obtenidos con CEINCI-LAB. Por esta razon no tiene 
mayor relevancia presentar el espectro inelastico de la Presa Minas-La Union-Pitahuina como 
tampoco lo tiene discutir el valor de reduccion de las fuerzas sfsmicas adoptado pero vale la 
pena indicar que un valor de R — 2 para una Presa de Tierra es alto. 
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5.4 Resultados del analisis sismico 

En las figuras 18 y 19 se presentan los desplazamientos y fuerzas estaticas 
equivalentes, obtenidas en los nodos, del analisis sismico. 



Figura 18 Desplazamientos modales inelasticos, en sentido horizontal. 



Figura 19 Fuerzas sismicas en sentido horizontal. 


5.5 Fuerzas Hidrostaticas 


Para las fuerzas hidrostaticas se utilizo la subrutina FUHIESTATICAS, la cual calcula 
una fuerza resultante perpendicular al piano inclinado en cada nodo de la cara aguas arriba. 
Luego saca las componentes horizontal y vertical. En la figura 20 se muestran los resultados en 
los que se han etiquetado los primeros valores correspondientes a la primera fila. 




X-HIDROSTATICAS 



Areas coloreadas de HIDROSTATICAS, X-HIDROSTATICAS. 


! 


34.737 

30.878 

27.018 

23.158 

19.298 

15.439 


11.579 

1 7.7194 
3.8597 
0 


Figura 20 Fuerzas Hidrostaticas en el sentido X. 


Estrictamente las fuerzas hidrostaticas se obtienen a nivel de la cara de la Presa pero el 
Programa FUHIESTATICAS reparte estas fuerzas en todos los nodos, considerando una 
variacion lineal de atenuacion de las fuerzas a lo largo de la Presa, es asi como en el borde 
extremo las fuerzas hidrostaticas son cero. 



Figura 21 Fuerzas Hidrostaticas en el sentido Y. 


En la figura 21 se presentan las fuerzas hidostaticas en sentido Y, actuando unicamente 
en el borde de la Presa mientras que las fuerzas en el sentido X se distribuyeron en los demas 
nodos. 
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Fuerzas Hidrodinamicas 


Las fuerzas hidrodinamicas fueron encontradas con el metodo de Housner. El metodo 
encuentra la presion hidrodinamica y las fuerzas horizontal y vertical en el borde de la cara en 
contacto con el agua. Para las fuerzas horizontales estas fuerzas se han distribuido en forma 
lineal a lo largo de la Presa, en forma similar a lo realizado con las Fuerzas Hidrostaticas. Los 
resultados se muestran en las figuras 22 y 23. 


X-HIDRODINAMICAS 



Areas coloreadas de HIDRODINAMICAS, X-HIDRODINAMICAS. 


Figura 22 Fuerzas Hidrodinamicas en el sentido X. 



Figura 23 Fuerzas Hidrodinamicas en el sentido Y. 


Fuerzas Totales 


Las fuerzas que intervienen en el sentido X son: las fuerzas sfsmicas, las hidrostaticas e 
Hidrodinamicas en X. En el sentido Y intervienen el peso propio, las fuerzas hidrostaticas y las 
fuerzas hidrodinamicas en esa direction. Estas fuerzas a nivel de nodos, son las que se 
presentan en las figuras 24 y 25 para el sentido X, Y, respectivamente. 



Figura 24 Fuerzas en sentido X., en los nodos. 



Figura 25 Fuerzas en sentido Y, en los nodos. 
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Desplazamientos Totales 

Con las fuerzas totales actuando en sentido X, Y, se hallan los desplazamientos totales 
a nivel de los nodos, previamente se encontro la matriz de rigidez asociada a todos los grados 
de libertad y se resolvio un sistema de ecuaciones lineales. Los resultados se muestran en las 
figuras 26 y 27, para el sentido X, Y, respectivamente. 



Figura 26 Desplazamientos totales en sentido X. 



Figura 27 Desplazamientos totales en sentido Y. 


6. COMENTARIOS Y CONCLUSIONES 

Si bien es cierto, el objetivo de este articulo era el mostrar como se utiliza el programa 
GID, en la etapa de post proceso, a partir de resultados obtenidos con los programas del 
sistema de computacion CEINCI-LAB , no es menos cierto que lo expuesto tiene el caracter 
general por lo que se debe proceder en forma similar con los resultados de otros programas. 

Se aspira que con los programas y explicacion presentada en este articulo el lector se 
encuentre en capacidad de presentar resultados elegantes, vistosos y claros, de problemas 
resueltos con elementos finitos, utilizando el programa GID como postprocesador. 
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